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Abstract 

We study the renormalization group evolution up to the fixed point of the lat- 
tice topological susceptibility in the 2-d 0(3) non-linear ci-model. We start 
with a discretization of the continuum topological charge by a local charge 
density, polynomial in the lattice fields. Among the different choices we pro- 
pose also a Symanzik-improved lattice topological charge. We check step by 
step in the renormalization group iteration the progressive dumping of quan- 
tum fiuctuations, which are responsible for the additive and multiplicative 
renormalizations of the lattice topological susceptibility with respect to the 
continuum definition. We find that already after three iterations these renor- 
malizations are negligible and an excellent approximation of the fixed point 
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is achieved. We also check by an explicit calculation that the assumption of 
slowly varying fields in iterating the renormalization group does not lead to 
a good approximation of the fixed point charge operator. 

PACS numbers: 05.50.+q 11. 15. Ha ll.15.Tk 75.10.Jm 
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I. INTRODUCTION 



In the field theoretical approach the lattice is an ultraviolet regulator for continuum 
quantum field theories, with the lattice spacing acting as cutoff; physical observables are 
constructed on the lattice in terms of local operators having the correct classical contin- 
uum limit. Two classes of effects hamper the extraction of physics from the lattice in this 
approach. First, discretized quantities are connected with the continuum ones through non- 
trivial renormalizations [I]; these renormalizations should be accurately evaluated. Second, 
a slight cutoff dependence of lattice quantities, not subtracted by the renormalization pro- 
cedure and vanishing in the a — )• limit, induces systematic errors (the well known cutoff 
effects) in the determination of continuum quantities. Both classes of effects are artifacts of 
discretization. 

A very elegant and radical solution to these problems comes from renormalization group 
(RG) theory [2]: all lattice actions and operators corresponding to the renormalized trajec- 
tory (i.e. the asymptotic fiow line in the space of couplings under repeated RG transforma- 
tions) are perfect, in the sense that they are free from lattice artifacts. Using perfect actions 
and perfect operators in numerical simulations, information about the continuum can be 
directly extracted from the lattice, since renormalizations and cutoff effects are absent. The 
first step in this ambitious program consists in constructing, for asymptotically free theories, 
fixed point (FP) actions and operators, which are perfect in the classical limit — > (being 
g the asymptotically free coupling). A formal argument of perturbative RG implies that FP 
actions are I-loop (quantum) perfect [3,4]. 

A method for the determination of FP actions for asymptotically free theories has been 
proposed in a recent paper [5]. The procedure has been applied to the 2-d 0(3) non-linear 
cr-model on the lattice and a parametrization of the FP action suitable for numerical 
simulations has been found. In Ref. [6] a numerical implementation of RG has led to the 
definition of a FP topological charge operator in the geometrical approach [7] for the same 
model. 
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In this paper we build a FP topological charge operator for the 0(3) cr-model in the 
framework of the field theoretical method [1]; this operator is used to extract the topological 
susceptibility of the model through numerical simulations. We adopt the FP parametrized 
action proposed in Ref. [5]. 

The paper is organized as follows: in Section II we review the techniques used in the 
construction of FP actions and operators and make some comments about the locality of 
renormalized and FP density operators; in Section III we outline the analytical determina- 
tion of the FP topological charge at lowest order in a slowly varying field approximation; 
in Section IV we describe the numerical procedure which implements the RG transforma- 
tions; by iterating this procedure we construct the RG fiow of topological charge operators 
starting from a field theoretical definition and eventually converging to the FP operator; in 
Section V we outline the prescriptions of the field theoretical method in the computation of 
the topological susceptibility; we determine numerically its multiplicative renormalization 
and perturbative tail, when the first three renormalized charge operators of the RG fiow are 
used; in Section VI, we study the scaling of the topological susceptibility by implementing 
the RG formalism in the usual MC techniques; in Section VII, we draw some conclusions. 



In order to clarify the further exposition and fix the notation, we briefiy review the RG 
techniques used in constructing FP actions and operators for asymptotically free theories. 
Let us consider an asymptotically free lattice spin theory. We write its partition func- 



where ^ is a spin field living on a lattice with spacing a, and is some lattice regular- 

ization of the continuum action. In the case of the 0(3) non-linear cr-model [dcf)] would be 
the 0(3)-invariant measure 



II. THE FP TOPOLOGICAL CHARGE OPERATOR 



tion: 
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n 

and a regularization of the continuum action 

/3Aont = ^/ (Px d^(p{x) ■ d^(p{x) , = 1 . (3^ 

A RG transformation with scale factor A for the spin system can be defined in the following 
way: 

f^'A'm _ l\in f3{Al<t,]+Tl<S>,<t,]) . 



here $ is the blocked spin living on the lattice with spacing a' = A a, related to a local 
average of the original spin variables; T[$, (p] is the blocking kernel, normalized in order to 
keep the partition function invariant under the transformation. The transformation used in 
Ref. [5] for the 2-d 0(3) cr-model is obtained by dividing the original square lattice into 
2x2 blocks and assigning to each of them a new block spin $„g . This spin is constructed 
by averaging the four original spins of the block. In the /3 — )• 00 limit, Eq. (4) can be solved 
in the saddle point approximation^: 

A'm=rnin{A[<P] + T[<^,<P]} . (5) 

{<p} 

In this limit the blocking kernel T[$, (p] turns out to be [5]: 

T[$,<^] = -«^[$„, <Pn\] , (6) 

n-g ragriB nEn-Q 

where {n : n G ^b} are the sites n of the original lattice belonging to the block labeled 
with ub] /€ is a free parameter. The fixed point A* of Eq. (5) is the starting point, on the 
/3 = 00 critical surface, of the renormalized trajectory, which defines lattice actions free of 
cutoff effects. As a consequence, the EP action = is perfect in the classical limit 

/3 — )• 00, i.e. it is the classical perfect action. 



1/3' = (3 — 0(1) owing to asymptotic freedom. 
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For a given configuration {$} the value of tlie renormalized action ^'[$] can be deter- 
mined by finding tlie configuration {^}, defined on tlie finer lattice, which minimizes the 
right hand side of Eq. (5). Furthermore, if the RG transformation is iterated k times, the 
equation corresponding to Fq. (5) is 

= min {^(°)[#)]+T[#-^),#)] + ... + T[$,(^(^)]} , (7) 

where ^('^^ is the field at the A;-th fine level; starting on the finest lattice with any regular- 
ization of the lattice action, a good approximation of the FP action can be obtained 
even for small values of k [5] . 

In the case of the 0(3) cr-model the parametrization used in Ref. [5] for is: 

ApM = /3 - <l>n ■ <l>n+r) (8) 

I ^ n,r 

rii ' 4'n2 

ni ,n2 ,ns,rn 

where p represents the perfect discretization of the Laplacian; the coefficients c(ni, ^2, ^3, ^4) 
and p[r) can be determined analytically by expanding the FP equation for slowly varying 
fields [5]. In the following we will always refer to the 24-couplings parametrization given 
in Table 4 of Ref. [5], where the couplings above the quartic order in the fields have been 
determined numerically by a fit procedure on the functional ^*[^], implementing Fq. (7) for 
an ensemble of configurations {$} . 

The RG transformation for an operator 0[(p] is defined as follows: 

0'[$] e-'^'^'I^] = J [d^] 0[^] e-/^(^M+^[^'^]) ; (9) 

in the /3 — )• 00 limit it is straightforward to derive: 

o'm = o[<i>^^[^] , (10) 

where { ^min[$] } is the solution of the saddle point equation (5). Fq. (10) says that the 
renormalized operator O' calculated on the blocked configuration {$} has the same value as 
the original operator O calculated on the minimizing configuration { ^min[$] }• The iteration 
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of Eq. (10) yields a succession of operators O^'^^ starting from some lattice regularization 
0^°^ satisfying the recursive equation: 



(11) 



The solution of (11) is: 



(12) 



being {^('^)[$]} the minimizing configuration on the finest lattice in Eq. (7). The fixed point 
operator Off is the k ^ oo limit of O^'^K A numerical determination of this limit can be 



practically viable if the EP operator is approximated by the renormalized operator after 
few iterations. Using the same RG argument of the EP action, one concludes that the EP 
operator is the classical perfect operator and its classical properties are the same of the 
corresponding continuum operator. 

In the following we will be concerned with the topological charge operator of the 0(3) 
non-linear cr-model; however the discussion is quite general and can be extended to every 
operator which can be expressed in terms of a local density. 

Let us suppose to be on the fixed point A* of the RG transformation (5). The minimizing 
configuration {^min[$]} on the fine lattice can be parametrized [6] in the following way: 



(^min)„ = X1"(^'^b) <i>nB + /3{n,lB,rnB,nB){l - ■ ^nB)^lB + ■ ■ ■ ■ (l^) 



Since we sit on the fixed point, the form of the parametrization is independent of the RG 
step. We define J(n, me) as the distance in coarse lattice units between the site niB and the 
site ub to which n belongs. We define the spread of the parametrization (13) as the maximum 
distance J(n,mB), over all the me's in the coarse lattice involved in the parametrization of 



discretized local density operator C'[^](n) as the maximum distance J(n, n') over the sites n' 
involved in the definition of C'[^](n). By definition of local operator, this spread must keep 
finite in lattice units as the continuum limit is reached. 
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Let us now consider Eq. (11) for the topological charge and suppose that Q^^ is 
expressed in terms of a local density operator: 

= . (14) 

n 

We can rewrite the last equation performing a blocking on the last summation: 

Q^'^[$] = E(E Q^'-'\<t>^Amn)]=Y.Q^'\mnB). (15) 

n-B VragriB / n-g 

The last equation defines implicitly the density Q^^\nB)- Suppose that the field parametriza- 
tion has spread r (in units of the coarse lattice spacing) and that Q^^~^\n) has spread 
Sk-i (in terms of the finer lattice spacing). Then, it is quite clear that Q^^^u-q) has 
spread ~ r + From this follows that, starting from a lattice regularization 

= En with spread s, every element of the succession Q^^'> can be parametrized 

in terms of a discretized local density having spread 

Sk ^ r + - + + + -r . (16) 
" 2 2 2 

Since ^min[$] can be parametrized in a local way [6], the operators Q^^\ A; = 1,2, . . ., and 
ultimately the FP charge Qfp itself, will still be expressed in terms of local densities Q^^\n) 
with finite spread, the limiting value b emg -Spp ~ 2r. This feature allows the application of 
heating techniques to determine the behavior of the renormalizations of the lattice operators 
Q^*^) (see the following). What however is not under control is the number of powers of the 
fields contained in the density operators. Indeed, due to the non-linear terms in the field 
parametrization (13), the iteration of Fq. (11) is likely to give rise to an uncontrolled growth 
of the powers of the fields involved in the parametrization of Q^^\n) as k grows. That Qyv 
in particular cannot be well approximated by a polynomial expression in the fields, can be 
concluded using the following two arguments. 

i) The FP charge operator has, as in the continuum, a discrete spectrum. This is clear 
considering Fq. (12) in the case O = Q: in the limit k ^ oo the minimizing configuration 
on the finest lattice ^('^)[$] approaches a classical continuum solution [6], on which every 
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starting regularization of the charge Q^^^ gives an integer value (the same for all). This 
property is consistent with the RG argument which ensures that Qpp is a perfect classical 
charge operator: the FP charge reflects, up to a minimum size, the discrete classical spectrum 
of the continuum operator 

Q = I (fx Q{x) , Q{x) = -!-e^j. e,jkMx)9,i<Pj{x)d^<pk{x) ■ (17) 

ii) On the lattice, differently from the continuum, each configuration can be continuously 
deformed into any other^. Joining arguments i) and ii) one concludes that Qfp cannot 
have an analytical functional dependence in the lattice field"^: one can produce for instance 
a discontinuous jump of Qpp by continuously deforming a Qpp = configuration into a 
Qfp = 1 one. 

III. Qfp for slowly varying fields 

A possible way to construct FP actions and operators, is to solve Fq. (11) making an ex- 
pansion for slowly varying fields. The parameter of this expansion is the mean angle between 
the spins in the configuration {(f>}. In this Section we present the analytical determination 
of (5fp[^] at lowest order. A necessary step in the calculation is finding the minimizing 
configuration { ^min[$] } at lowest order in the expansion, i.e. the term linear in $ in the 
parametrization (13). This amounts to solving the saddle point equation (5) in the free field 
theory case [8,5]. We express a(n,nB) of Fq. (13) in Fourier transform: 

a(n,nB) = ^ e''^'^"'^^^'' ' (18) 

where the subscript B means that the integration domain is the first Brillouin zone 
[— 7r/a, TT / a] X [— vr/a, vr/a] and both n and ub are expressed in the coordinates of the coarse 

^On the continuum this is impossible for starting and final configurations having different topo- 
logical charge. 

■^In this respect, Qpp resembles the geometrical definition of the lattice topological charge [7]. 
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lattice. The result can be written as follows: 

F{k) = -a{k)a'{2k)-^ H 2 cos ^ ; (19) 

here, a[k) is the inverse of the Laplacian p[k) and the prime means after one RG transfor- 
mation. 

We used the following 0(3)-invariant odd-parity parametrization of the topological 
charge operator at lowest order in the fields: 

SttQ = /"(^i -^2,ni - ns) e,jk M'^i) ^^(^2) ^^(^^3) (20) 

ni,n2,n3 

JB (27rj (27rj 
^) 2. L ^2 ^2 M(gi + mi7r,g2 + m27r) 

^ ^ mi,m2 "^^ " " 

X e^jfc (|)^{-ql - rriiTT - q2 - m27r) <pj{qi + mivr) ^^(52 + ?7^27r) , 

where mi, m2 are 2-vectors whose components take the values and 1; the second equation 
follows from expressing the summation over the sites n as a summation over the blocked sites 
ub and over the four sites of the block defining the fixed ne; B' means the "halved" Brillouin 
zone [— 7r/(2a), 7r/(2a)] X [— 7r/(2a), 7r/(2a)]. If we now write Q after one RG transformation 
as 

sttQ' = [ f^ff i,'{q,^q,)e,,kH-qi-q2)'^Aqi)'^k{q2) (21) 

and, according to Eq. (14), plug the lowest order approximation of ^min[$] (in Fourier 
transform) in the last integral in Eq. (20), we find: 

/(2gi,2g2) = 229I43 m(9i + "ii7r,52 + m27r) (22) 

mi,m2 

p'{2q,+2q2)p'{2q,)p'{2q2) 



p{qi + miTT + q2 + m2Tr)p{qi + miTr)p{q2 + m27rj 



nf qi + niiTT + q2 + m2Tr\ -pf ^ f qi + miTr\ ^ f q2 + m2Tr 
2cos [ 2 j^ll^cos 

The EP //-function, 

/^FP 1 C9jI1 be found iterating the previous relation; it comes out that 
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Mfp (91,92)= /«cont(9i + 2mi7r,52 + 2m27r) (23) 

mi ,m2 €-2"^ 

/Ofp(9i + 92)/Ofp(9i)/Ofp(92) 



X 



[qi + 2mi7r + 92 + 2m27ry{qi + 2mi7ry{q2 + 2m27r)2 

sin((5i + 52)/2)^ -pr sin(5i/2)j, ^ sin(52/2)a 



V" ( (?i + ?2)/2 + miTT + m27r)^ V" ( 9i/2 + niiTr)^ V" ( ^2/2 + m27r)<, 
where /Ucont(p, <?) = ^iivPii(lv and, using /Jfp, we have assumed that the action is also at the 
FP. We give for completeness the expression of p^-p [8,5]: 

-1 = V ^ TT siii^(PA^/2) ^ . . 

k^{p+m' ■^/(p./2 + g^ ^ 3« ■ ^ ^ 

As in [5] we use n = 2. In Table I we show the coefficients //(m, n) up to 0(10"^) (coefficients 
which can be related to these by symmetry transformations are omitted). The property of 
locality of //(m, n) is consistent with the general discussion of previous Section about the 
spread of FP operators. 

We have tested the leading Qfp operator in the slowly varying field expansion, observing 
that it fails to reproduce the integer-valued spectrum of the continuum even on quite smooth 
configurations. We conclude that even on these configurations the neglected higher order 
operators dominate the expansion. This agrees with the theoretical conclusion at the end of 
the previous Section. 



IV. NUMERICAL DETERMINATION OF Qyv 

An alternative way to face the problem of the determination of the FP charge of a given 
lattice configuration, is to solve the FP equation for the charge operator numerically. This 
approach has been described and applied recently to the 0(3) cr-model in Ref. [6] for the 
geometrical definition of the lattice topological charge. Starting from a certain configuration 
{$} defined on the "coarse lattice" L^^\ the minimizing configuration { ^('^)[$] } of Fq. (12) 
is singled out by constructing the succession of configurations defined on finer and finer 
lattices: 

{$}, ... , {#)}^ {#)[$]} ; (25) 
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given the configuration {^('^ on the lattice L^'^ ^\ the configuration { ^('^^ } on the lattice 
iv^*^) is determined by numerically minimizing the quantity 

A*[<P']+T[<P^''-'\<P'] (26) 

over all the configurations {(f)'} living on the lattice L^'^K After few steps of this iteration, 
the irrelevant content in = ^^'^^[$] ] (Q^^^ is the starting lattice operator) is 

negligible, and this quantity is a good approximation of (5fp[$]- We have performed this 
program using three different starting lattice regularization Q^^^ of the charge: 
(1) the naive regularization Qst^'- 

1 



Stt 



St v^) = ^ e^i' e^jfc (f)^{x)V f,(f)j{x)V ^(f)k{x) , (27) 



where is the symmetrized lattice derivative 



Vi^(l),{x) = ^(^j(x + /w) - (t)tix - 1^)) ; (28) 

(0) . 



(2) the Symanzik tree-level improved (up to 0[a )) operator Qsym- 



)(0) 
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where is the improved lattice derivative: 

1225 245 49 5 , , 

Da = V V + V V J) , 30 

^ 1024 " 1024 " 1024 " 1024 " ^ ^ 

V|,''V«(^) = —i(t)tix + nii)-(l),{x-nii)) ; 
Zn 

(3) the FP operator at lowest order in the slowly varying field approximation Qgyp derived 
in the previous Section'*. 

We performed the minimization of the quantity in Eq. (26) by a local Metropolis algo- 
rithm accepting only lowering changes. The algorithm started from the minimizing config- 
uration at lowest (linear) order in the slowly varying fields expansion: 



Cf,^ e,jk Mx)D,,<Pj{x)D^<pk{x) , (29) 



"^Using the first 12 leading couplings. 
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(31) 



KB 

The number of minimization sweeps needed to reach the minimum decreases very rapidly 
with k and at the third step no further minimization beyond the hnear approximation is 



The starting topological charges give clearly different results at the step zero, since large 
cutoff effects are present which are unequally handled by the various definitions; but, under 
iteration of the RG transformation, these UV effects are gradually erased and the three 
different field theoretical definitions of the lattice charge operator converge to the same 
operator Qpp. 



The problem of extracting from the lattice the topological susceptibility of the 0(3) cr- 
model has been studied for many years [6,7,9-15] using different approaches. The topological 
susceptibility has the dimension of a mass squared, so admitting that it is a physical quantity, 
i.e. renormalization group invariant, its (adimensional) lattice regularization should scale 
in the continuum limit with i^"^, where ^ is the correlation length in lattice units. In fact, 
this scenario is obscured by a prevision of the semiclassical approximation [16] and some 
numerical evidences [17,18] which indicate that the topology of the model is UV dominated. 

In the field theoretical method [1] a lattice topological charge density operator is de- 
fined as a local operator Q{x) having the appropriate classical continuum limit. The matrix 
elements of Q{x) are related to those in the continuum through a finite multiplicative renor- 
malization: 



The lattice version of the topological susceptibility is simply given by x = 1 / L'^ (( Q{^) Y) 
[L is the lattice size). It is connected to the continuum counterpart by the relation 




V. THE FIELD THEORETICAL METHOD 
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x(/3) = a'^(/3)'xcont + + 0{a^) ; (33) 

where P{fi) is the perturbative expectation value of the lattice topological susceptibility 
(the so called perturbative tail)^. Both Z{fi) and P{fi) are artifacts of discretization in- 
duced by the lattice quantum fluctuations. Eq. (33) allows to extract Xcont, the topological 
susceptibility of the continuum theory, from MC determinations of the lattice regularized 
susceptibility x- An intermediate step is the evaluation of the renormalizations Z{fi) and 
P{fi)- In standard perturbation theory they are calculated as power series in [13]; an 
alternative approach is the Monte Carlo "heating" technique [19,13,18] (see the following). 

Let's now come to the central issue of this paper, i.e. the application of the field theo- 
retical method in combination with the FP lattice regularization of the topological charge. 
We move towards the FP operator through the operators Q^^"^ [k = 1,2,3) and construct 
the lattice susceptibilities x^*^^ = (^{ J2xQ^'^H^) )^)- The operators Q^*^^ still satisfy (as 

the starting operator Q^^^) the prescriptions of field theory, since they can be expressed in 
terms of a local charge density with finite spread (see Section II). This allows us to write: 

xW(/3) = a2zW(/3)2 Xcont + P^'Hf){I) + 0{a') . (34) 

The evaluation in perturbation theory of the renormalizations Z^'^'^ and P^*^) is not practicable 
in this case because the analytical dependence of the operators Q^*^) on the field (f) is not 
known; consequently, we turn to the numerical approach: the heating method. 

The heating method [20,19] consists (for details see also [13,18]) in constructing on the 
lattice ensembles of configurations {Cj}, each configuration of the ensemble being obtained 
by performing a sequence of t local Monte Carlo sweeps starting from a discretized classical 
configuration Co - a large instanton or a fiat configuration. For small t, the heating process 
thermalizes only the small-range fiuctuations which are responsible for the renormalizations: 



^In the r.h.s. of Eq. (33) the mixing of the topological susceptibility with the action density has 
not been considered, since it is numerically negligible [18,15]. 
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when the starting configuration is a large instanton (fiat configuration), measuring Q (or 
X) on the ensembles {Cj} a plateau at the value of Z{fi) {P{fi)) is observed after a certain 
time, not depending on /3, corresponding to the thermalization of quantum fiuctuations. The 
adaptation of this technique to the present case needs the building of the ensembles {C^*^^} 
at the three fine levels. Starting from a configuration on the coarse lattice {$} G {C**}, 
the corresponding configuration at the k-ih fine level {^('^)[$]} G {Cj*^^} was obtained by 
performing the minimization procedure described in Section IV. Averaging the starting 
definition of the topological charge (susceptibility) on the ensembles {C^*^^} amounts to 
performing the heating procedure for Q^^"^ (x('^)),so that the renormalization constant Z^^"^ 
{P(^y) at each RG level can be determined. In Tables II - III we show the values of Z^'^'^ and 
P^^^ k = 1,2,3 (in the case /3 = 1.15 also k = 0) for five different /3's with the Symanzik 
improved topological charge as the starting operator: it is interesting to observe that Z and 
P converge rapidly towards 1 and respectively when k increases, thus indicating that the 
FP topological charge and susceptibility are insensitive to quantum fiuctuations. 

In Figs. 1 and 2 we show the behavior of ((5^°^)(/3 = 1.15) and {Q^^'^){/3 = 1.15) during the 
heating of an extended 1-instanton (with fixed boundary conditions). The most important 
observation is that the time of thermalization of the quantum fiuctuations of Q^^^ - not yet 
erased by the RG procedure - is very short (~ 4 heating steps), thus indicating that the 
spread of this operator is small (of the order of few lattice spacings). This gives a numerical 
support to the arguments of Section II about the spread of renormalized operators defined 
starting from a local density. 

From the analysis of the data of heating we conclude that Qg'y^^ is the starting operator 
closest to the fixed point, since it yields at a given k the values of Z and P closest to 1 and 
respectively (see Table IV for Z(/3 = 1.15)) ^. In the following we will be concerned about 



^That Q 

Sym is closest to the FP operator could support the idea that Symanzik's improvement 
program goes in the same direction of the RG improvement. 
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the data coming from this starting operator. 



VI. MC SIMULATIONS 

We performed Monte Carlo simulations on a starting "coarse" lattice 90 X 90 for five 
values of /3 (1., 1.05, 1.10, 1.15, 1.20), ^ ranging from 12 to 35 lattice units. We measured 
Qsym (oiily ill the case [3 = 1.15), Qsym? Qsym ^iid Qsym usiug the numerical technique 
discussed in Section IV. We considered ~ 5000 decorrelated starting configurations at 
thermal equilibrium for each value of /3. 

In Table V and in Fig. 3 we present values of Xcont ■ Cp? where = i^a, at the five 
/3's considered on each RG level: Xcont has been extracted according to Eq. (34); ^ was 
determined by evaluating the exponential decay of the wall-wall correlation function on 
lattices satisfying the condition Z/i^ ~ 7 (in order to keep safe from finite size effects)''. 

Data in Table V allows two observations: the vertical reading shows that there is no 
scaling up to correlation length ^ ~ 35; the horizontal reading shows that the physical 
signal Xcont increases after each RG step. The first observation is a consequence of the 
UV dominated nature of the topology of the model, already well established in literature; 
the second effect is subtler, related to the different scale invariance properties of the three 
operators Q^'^^ k = 1,2,3. Being classically perfect, Qfp should have a scale invariant 
spectrum, so it should attribute the correct topological charge to lattice configurations up 
to a very small size, close to the lattice unit. In the case of Q^*^) this property is only 
approximate, improving for increasing k. This is visually clear in Fig. 4, where we show the 
behavior of topological charge at the coarse, the first fine and the second fine levels on small 
size classical configurations (1-instantons with fixed boundary conditions). We can see that 
the topological charge at the coarse level underestimates the correct continuum value Q = 1 



''We are indebted to the authors of Ref. [6] for allowing us to use their cluster algorithm program 
for measurements of correlation length. 
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on classical configurations with small size, while the correct value is asymptotically attained 
as k increases. We conclude that the slight dependence of our determination of Xcont on the 
"fineness level" of the lattice, A;, is an effect of the progressive saturation of the topological 
signal coming from small size instantons. We observe that at A; = 3 this saturation can 
be considered practically complete; we argue therefore that only a negligible improvement 
would come from further RG steps. It is also evident from Fig. 4 that no topological signal 
is observable on the lattice below 0.9 lattice units, even when using a FP operator and a 
FP action. Since most of the topological signal in this model comes from the region of 
small lengths, a cutoff effect is still present in the lattice theory, which manifests itself in 
the unphysical dependence of Xcont on ^ (non-scaling). 

VII. SUMMARY AND CONCLUSIONS 

The theory of Wilson's renormalization group states the existence of perfect lattice oper- 
ators, i.e. operators discretized on the lattice which are free of cutoff effects. These operators 
correspond to the renormalized trajectory of a given renormalization group transformation. 
The fixed point of the transformation defines the fixed point operator, the classically perfect 
lattice operator. 

In this paper we have studied the evolution towards the fixed point of the two-point 
correlation function of the topological charge, i.e. the topological susceptibility, in the 2-d 
0(3) non-linear cr-model on the lattice. For this purpose, we have discretized the continuum 
topological charge on the lattice in terms of a local operator, polynomial in the lattice fields, 
and we have implemented the renormalization group procedure in a numerical fashion. 

We have observed that the quantum fiuctuations are progressively erased as the fixed 
point is approached, thus leading to an integer-valued spectrum of the topological charge 
on the lattice. The amount of quantum fiuctuations is quantified by the renormalization 
constants Z[l3) and P{/3) entering the lattice definition of the topological susceptibility. We 
have measured these constants by the heating method observing their rapid disappearance 
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with the iteration of the renormahzation group transformation. The operators obtained 
at each step of the renormahzation group procedure can be considered as improvements 
of the starting operators, in the sense of a suppression of the quantum noise around the 
physical signaL A similar improvement, consisting in the reduction of the renormahzation 
constants, was obtained in Ref. [21] by using smearing techniques. We have also calculated 
the physical value of the topological susceptibility at thermal equilibrium by subtracting the 
effect of the renormalizations from the lattice signal at each renormahzation group step. The 
physical quantity obtained increases with the renormahzation group step up to convergence 
within a few percent of accuracy. This convergence corresponds to the saturation of the 
topological signal by smaller and smaller instantons, down to the minimum size allowed on 
the lattice. The physical signal at the hxed point does not exhibit a scaling behavior, since 
the topological contribution of instantons with size lower than the critical size, which is still 
lost on the lattice, is dominating. 

A procedure similar to the one described in this paper has already been applied in 
Ref. [6] using the geometrical approach to the lattice topological charge: we observe that the 
two different methods, while diverging at the starting point, give asymptotically consistent 
results when the renormahzation group procedure is carried on. 

We have also presented an analytic determination of the hxed point topological charge 
in the approximation of slowly varying helds, by iterating the renormahzation group trans- 
formation at the leading order. The resulting operator fails to reproduce the integer-valued 
spectrum of the continuum even on quite smooth conhgurations, thus leading to the conclu- 
sion that the neglected higher order terms are important. The negative result obtained is 
not completely unexpected since an analytic expression in the lattice helds for the topologi- 
cal charge operator cannot reproduce a discrete (and, therefore, discontinuous) spectrum on 
the lattice, where each conhguration can be continuously transformed into any other. 

The present paper can be considered as a preliminary work in view of the determination 
of the hxed point topological charge for the SU(3) gauge theory. For this theory, however, 
the numerical procedure applied in the present paper is no longer viable, owing to the huge 
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computation time needed for the determination of the configurations on the finer lattices. 
We retain that a promising approach is to search for a (few-couphngs) parametrization of 
the fixed point topological charge operator which works even on coarse lattices. 
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FIGURES 

FIG. 1. Topological charge at the coarse level during the heating of an extended 1-instanton 
(with fixed boundary conditions) at /3 = 1.15. The line represent the fit (with errors) of the data 
on the plateau. 

FIG. 2. Topological charge at the first fine level during the heating of an extended 1-instanton 
(with fixed boundary conditions) at /3 = 1.15. The line represent the fit (with errors) of the data 
on the plateau. 

FIG. 3. X ' &p i^p = ^(i) foi" five correlation lengths ^ at the various RG steps. The starting 
topological charge regularization is the Symanzik charge. Data at the coarse, 2nd and 3rd fine 
level have been slightly shifted for the sake of readability. 

FIG. 4. Topological charge at the coarse, first fine and second fine levels on classical configu- 
rations consisting of 1-instantons (with fixed boundary conditions) with different size. Instantons 
have been smoothed out by some cooling steps before the charge measurement. 
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TABLES 



TABLE L Largest non-equivalent values of ii{m,n). 


m 


n 




li{m, n) 


m 


n 


li{m, n) 


(1,1) 


(1,0) 




5.00100x10-2 


(3,0) 


(2,1) 


-4.35376x10-^ 


(2,0) 


(1,1) 




-6.27316x10-3 


(3,1) 


(2,0) 


4.15550x10-5 


(2,1) 


(1,1) 




-4.90036x10-3 


(3,2) 


(2,2) 


3.90132x10-5 


(2,1) 


(2,0) 




1.52221x10-3 


(3,2) 


(2,1) 


-2.53691x10-5 


(2,1) 


(1,2) 




-9.39105x10-'* 


(3,1) 


(2,2) 


2.20664x10-5 


(2,1) 


(2,-1) 




3.55142x10-'* 


(2,2) 


(2,1) 


1.51901x10-5 


(2,2) 


(2,0) 




1.16270x10-'* 


(3,1) 


(3,0) 


1.37893x10-5 


(3,1) 


(2,1) 




-8.39146x10-^ 


(3,2) 


(3,1) 


-1 981 19x1 0-5 


TABLE IL Z{p) for the different /3's at the various 


RG levels, k. The starting 


operator is Q^^^. 






coarse 


1. level 




2. level 


3. level 


1.00 
1.05 
1.10 
1.15 
1.20 




0.450(6) 


0.892(9) 
0.910(8) 
0.921(6) 
0.929(6) 
0.936(5) 




0.980(5) 
0.990(3) 
0.9929(23) 
0.9931(21) 
0.9939(15) 


0.9950(27) 
0.9999(11) 
0.9992(8) 
0.9990(6) 
0.9991(3) 
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TABLE III. P(/3) X 10'* for the different /3's at tlie various RG levels, k. The starting operator 



(3 coarse 1. level 2. level 3. level 

1.00 0.571(22) 0.177(7) 0.0466(24) 

1.05 0.179(10) 0.0276(22) 0.0037(5) 

1.10 0.123(7) 0.0241(19) 0.0022(3) 

1.15 0.804(13) 0.094(5) 0.0116(10) 0.00064(12) 

1.20 0.066(3) 0.0088(9) 0.00061(11) 



TABLE IV. Z(P) at /3 = 1.15 for the three different starting regularizations of the topological 
charge. 

coarse 1. level 2. level 3. level 

git^ 0.411(4) 0.855(5) 0.9625(25) 0.9903(10) 

gg^^ 0.450(6) 0.929(6) 0.9931(21) 0.9990(6) 

g^°^p 0.384(3) 0.815(5) 0.9478(27) 0.9860(11) 
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TABLE V. Topological susceptibility in physical units up to the third RG level (^p = ^a). For 
comparison, we also report in the last column the data of Ref. [6]. 



P S Xcoarse ' Xl. level ' X2. level ' X3. level ' Sp Xl. level ' \- \ 

1.00 12.16(3) 0.1448(15) 

12.27(3) 0.137(7) 0.150(6) 0.151(5) 

1.05 15.81(6) 0.161(8) 0.172(6) 0.173(5) 

1.10 20.34(9) 0.173(8) 0.185(7) 0.188(6) 

20.40(9) 0.1893(27) 

1.15 26.26(14) 0.160(13) 0.177(8) 0.192(7) 0.195(7) 

1.20 34.13(25) 0.194(11) 0.210(10) 0.213(9) 

34.4(3) 0.224(5) 
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